Metastable Random Field Ising model with exchange enhancement: a simple model 

for Exchange Bias 
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We present a simple model that allows hysteresis loops with exchange bias to be reproduced. The 
model is a modification of the T — random field Ising model driven by an external field and with 
synchronous local relaxation dynamics. The main novelty of the model is that a certain fraction / 
of the exchange constants between neighbouring spins is enhanced to a very large value Je- The 
model allows the dependence of the exchange bias and other properties of the hysteresis loops to 
be analyzed as a function of the parameters of the model: the fraction / of enhanced bonds, the 
amount of the enhancement Je and the amount of disorder which is controlled by the width a of 
the Gaussian distribution of the random fields. 

PACS numbers: 75.40.Mg, 75.50.Lk, 05.50.+q, 75.10.Nr, 75.60.Ej 



I. INTRODUCTION 

Hysteresis and metastability are intriguing phenom- 
ena with implications in both fundamental and applied 
physics. Magnetic systems are the prototypical example 
of thermodynamic systems exhibiting hysteresis cycles 
for which different theoretical approaches have been pro- 
posed Q. Besides the classical micro-magnetic analysis, 
based on a continuous description of the magnetic prop- 
erties of the system, more recently a lot of efforts have 
been devoted to the study of lattice models. For example, 
the zero-temperature Random Field Ising model (RFIM) 
driven by an external field with convenient metastable 
dynamics has been very successful in qualitatively ex- 
plaining some basic properties of rate independent hys- 
teresis loops ||, |[ [| |j| . The most important achievement 
of this model has been to give a simultaneous explanation 
of the effect of disorder on the hysteresis loops and the 
existence of Barkhaussen noise with critical properties. 
Less attention has been paid to the use of such models for 
understanding other interesting features of the hysteresis 
loops such as remanence, coercivity, minor loop proper- 
ties or exchange bias (EB) || (?J, which is the property 
we will focus our attention on here. 

We present a modification of the zero-temperature 
RFIM that allows magnetic hysteresis loops with EB to 
be reproduced. The main characteristic of EB is that the 
hysteresis loops, represented as magnetization m versus 
external applied field H, are not centred on H = but ex- 
hibit a displacement in the field axis by an amount Heb 
(exchange bias field) |6|, |?| . This property has received a 
lot of attention recently, since the possibility of finding 
systems with large EB has enormous technological inter- 
est H. Experimentally EB has been found in different 
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magnetic systems |(| f?J. The basic ingredient for EB is 
the existence of interfaces between ferromagnetic (FM) 
and antiferromagnetic (AFM) systems, where coupling 
can be induced after field cooling from above the Neel 
temperature of the AFM. This heat treatment freezes 
some of the magnetic moments at the interface which are 
supposed to be responsible for the occurrence of EB. The 
prototype is a FM/AFM bilayer, for instance Co/CoO 
f|, NiFe/NiMn §J, Fe/FeF 2 §J and Fe/MnF 2 jjj. 
This effect is also observed in granular systems formed 
by small particles with a ferromagnetic core covered by 
their native antiferromagnetic oxide uM . 



Different models have been proposed to understand 
EB. Although a basic qualitative explanation was given 
40 years ago jl4j a deep understanding of the phe- 
nomenon has not yet been achieved Jl5|, Different 
features remain unclear: the role played by the AFM 
thickness ^ |T7| , the formation of domain walls (l^, [l9| , 
whether the frozen spins belong to the FM or to the AFM 
layer [^0), etc... Especially intriguing is the fact that EB 
not only occurs, in uncompensated AFM layers which 
exhibit a net magnetization after being cooled, but also 
in compensated AFM layers with zero net interface mag- 
netization ElL 



The aim of the present paper is to introduce a very 
simple model with a new mechanism for the explanation 
of EB in totally compensated layers. The model is based 
on a lattice spin system with metastable dynamics for 
which some of the exchange interactions show a marked 
enhancement. In section II the Hamiltonian and the de- 
tailed mestastable dynamics are presented. In section 
III we show the results of the numerical simulations. In 
section |fv| we discuss the possible physical origin of the 
exchange enhancement. In section |\] we compare with 
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available experimental data and, finally, in section |VI 
summarize and conclude. 
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II. MODEL 

The model is intended to reproduce the properties of 
the ferromagnetic layer only, and regards the AFM part 
to be totally quenched and thus does not contribute to 
the net magnetization (compensated AFM layer) . Conse- 
quently, we consider the 2d RFIM on the square lattice, 
although a generalization to bulk ferromagnets or thin 
layers could easily be implemented. Note, however, that 
the AFM plays an indirect role, as will be discussed later. 

The mathematical formulation of the model is very 
similar to the RFIM on a square lattice with size N — 
L x L. On each lattice site we define a spin variable 
Si which takes values ±1. The Hamiltonian, in reduced 
units, reads: 



stable situation is given by all the spins Si = 1. We de- 
crease the external field H until Hi vanishes on a certain 
spin. The spin is then reversed keeping H constant. This 
reversal may destabilize some of the neighbouring spins 
which arc then reversed simultaneously. This is the be- 
ginning of an avalanche. The avalanches proceed until 
a new stable situation with all the spins Si aligned with 
their respective local fields Hi is reached. We can then 
continue decreasing the external field H. 

Most of the calculated properties are averaged over a 
large number (~ 10 3 ) of different realizations of disorder. 
Averages are indicated by the symbol (•). We will con- 
sider, without loss of generality, that J = 1. Therefore, 
from now on, magnetic fields and energies are given in 
units of J. 



n 



N 



rjSiSj 



J2hSi-Hj2Si 



(1) 



The first sum is the ferromagnetic exchange contribution 
that extends over nearest neighbour pairs ( Jij > 0). The 
second sum accounts for the interaction with quenched 
random fields hi, which stands for the disorder present 
in any ferromagnetic system. The last term is the inter- 
action with the external driving field H. The random 
fields hi are independent and distributed according to a 
Gaussian probability density: 



p(hi 



1 



27TO- 



(2) 



where a is the standard deviation of the random fields 
and controls the amount of disorder in the system. The 
novelty of the model is in the values of the exchange con- 
stants Jij, which are not equal for all spin pairs: we con- 
sider that Jij — J except for a fraction / of the bonds 
(selected at random) for which = Jg >> J. This 
fraction of bonds is supposed to contain the effect of the 
quenched antiferromagnetic layer. A physical reason for 
this local exchange enhancement will be discussed in sec- 
tion IV. Such a bond distribution can be mathematically 



expressed as: 

p(Jy) = (1 - f)S(Jij ~J) + fS(Jij ~ J E ) (3) 

We have focused our study in the region of small values of 
/ (/ < 0.06). The magnetization of the system is defined 
as m — J2i=i Si/N. For the analysis of the hysteresis 
loops we use the so called synchronous local relaxation 
dynamics. This is the standard dynamics used in previ- 
ous studies of the zero-temperature RFIM . Each spin 
Si flips according to the sign of its local field Hi given 
by: 



Hi — JijSj + H + hi 

3=1 



(4) 



where the first sum extends over the four neighbours of 
Si. We start with a value of H large enough so that the 



III. RESULTS 

In Fig. [l] we present an example of a hysteresis loop 
obtained with the numerical simulation of a system with 
L = 50, J E = 20, / = 0.03 and a = 1.65. The external 
field is cycled between H — ±2.7. As can be seen, the 
loop exhibits remarkable EB. 




FIG. 1: Example of a hysteresis loop exhibiting exchange 
bias obtained from a numerical simulation of an L — 50 sys- 
tem with Je = 20, / = 0.03 and a = 1.65. The external field 
has been swept from 2.7 to -2.7. H c i and H C 2 indicate the 
coercive fields of the decreasing and the increasing branches 
respectively. 



At first glance it may look surprising that the model 
defined in the previous section displays such asymme- 
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try, since the Hamiltonian is totally symmetric under 
the changes Si — ► — Si and H — * —H. The reason is 
that the hysteresis loop shown in figure [l] corresponds, 
strictly speaking, to a minor loop. This is revealed in 
Fig. H where the system is cycled between H = ±19.4 
(which is a field that is one order of magnitude larger 
than the coercive field H c i). Due to the existence of a 
tiny fraction of very large exchange interactions, the to- 
tal loop exhibits long, flat plateaux in which the system 
behaves reversibly exactly as if it was saturated. Only 
when cycling between extremely large values of the exter- 
nal field, does one obtain the symmetric hysteresis loop. 
Therefore, the loops with EB are incomplete loops and 
are accompanied by a magnetization shift. 



m B 



m sl + m s2 1 + m s2 



(7) 




FIG. 2: The same example of Fig. 1 revealing the partial loop 
with exchange bias and the total loop between H = ±19.4 
which is symmetric. 



In order to perform a quantitative analysis of this be- 
haviour, averaging over different realizations of disorder 
is carried out. The hysteresis loops are systematically 
obtained according to the following protocol: decreasing 
the field from H = +oo to 2H c i < and increasing the 
field again to H — +oo, where H c \ is the coercive field 
in the decreasing branch. We also compute the pseudo 
coercive field H C 2 in the increasing branch (see figure [l]). 

The criteria for choosing the value 2iJ cl as a returning 
point is similar to the criteria used in many experimental 
cases. One could easily change this limit to ?>H c i or 4H c i 
without changing the results provided that Je is large 
enough. This can be easily understood from the flat tails 
in the full hysteresis loop shown in figure |[ 

The EB field Heb , the coercivity AH and the magne- 
tization shift mB are defined as: 



H 



cl 



H, 



c2 



Heb = 



AH = H C 2 — H c \ 



(5) 



(6) 



According to this definition, the loops shifted to the left 
on the H axis (as occurs with the loops in the present 
paper) will have negative exchange bias field. Figures || 
and [| show the dependence of (\H EB \) and AH on J E 
for a = 1 and different values of / as indicated by the 
legend. Even for very low values of /, (\Heb\) increases 
and saturates for large enough values of Je- 




FIG. 3: Dependence of the exchange bias (\Heb\) on the 
exchange enhancement Je for a — 1 and different values of 
/ as indicated by the legend. Data have been obtained by 
averaging 1000 realizations of a system of size L — 50. The 
lines are a guide to the eye. Statistical error bars are smaller 
than the symbol size. 



In the case of coercivity (Fig. ||) two important results 
should be emphasized: the increase (almost 40 % in cer- 
tain cases) in coercivity for intermediate values of Je and 
the saturation at a constant value (which depends on /) 
for large J E . Such limiting values at large J E however, 
are smaller than the coercivity of the system without ex- 
change enhancement. 



In order to analyze the dependence of the system prop- 
erties on the amount of disorder a, we choose a value 
of Je that is large enough so that (\H EB \) has reached 
the limiting maximum value (see Fig. ||). Figures |s| and 
^ show the behaviour of (\Heb\) and (AH) versus a 
for Je = 20 and different values of / as indicated by 
the legend. Unexpectedly, (\Heb\) shows non-monotonic 
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FIG. 4: Dependence of the coercivity (AH) on the exchange 
enhancement Je for a = 1 and different values of / as indi- 
cated by the legend. Data have been obtained by averaging 
1000 realizations of a system of size L = 50. The lines are a 
guide to the eye. Statistical error bars are smaller than the 
symbol size. 

behaviour with a, first decreasing until a minimum is 
reached, but which increases slowly for large amounts of 
disorder. 



The marked variation of (\Heb\) and (AH) for values 
of / between 0.015 and 0.025 when a — ► is associated 
with the fact that the ascending part of the hysteresis 
loop is very sensitive to the existence of nuclei of unre- 
versed spins in the negative magnetized state. For very 
low values of / we expect that all the nuclei of unreversed 
spins will be formed by two positive spins joined by an 
enhanced bond. The negative spins surrounding such a 
nucleus will flip (in the a — > limit) around H = 2 (see 
Eq. Q). However, for larger values of /, larger nuclei 
will exist. For instance, a nucleus formed by three spins 
joined by two perpendicular bonds, both with exchange 
enhancement, acts as a nucleating seed which triggers 
the avalanches towards the positive magnetization phase 
when H = 0. When / is large enough, such that the 
probability for such nuclei is significantly different from 
zero, the coercive field for the ascending branch decreases 
from 2 to 0, thus increasing Heb and decreasing AH. 



The non-monotonic behaviour of the EB when disorder 
is increased better seen by plotting (\Hb\) / (AH), which 




FIG. 5: Dependence of the exchange bias (\Heb\) on the 
amount of disorder a for Je = 20 and different values of 
/ as indicated by the legend. Data have been obtained by 
averaging 1000 realizations of a system of size L — 50. The 
lines are a guide to the eye. Statistical error bars are smaller 
than the symbol size. 

is a dimensionless quantity and is more relevant from the 
experimental point of view. Note that for large / and a, 
we can find strongly biased hysteresis loops for which H c \ 
and H C 2 are negative. An example, obtained by sweeping 
the field between H = ±4.5 is shown in figure ||. 



IV. DISCUSSION 

One of the reasons for the present work was to ob- 
tain hysteresis loops with EB by modifying the zero- 
temperature RFIM as little as possible. The straight- 
forward, naive idea would be to consider non-symmetric 
distributions of random fields, for which the average 
(hi) 7^ will create the displacement of the loop. From 
our point of view, this will correspond to the effect of 
an uncompensated AFM layer. Constraining ourselves 
to the inclusion of compensated disorder, we have found 
that any symmetric distribution of random fields can- 
not give EB. To understand this, suppose that a certain 
fraction of spins is pinned by very positive and (symmet- 
rically) very negative random fields. Fig. |](a) shows the 
schematic hysteresis loop corresponding to such a sys- 
tem. The spins with more positive random fields, which 
are the last to reverse in the decreasing branch, will be 
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FIG. 6: Dependence of the coercivity (AH) on the amount 
of disorder a for Je = 20 and different values of / as indicated 
by the legend. Data have been obtained by averaging 1000 
realizations of a system of size L = 50. The lines are a guide 
to the eye. Statistical error bars are smaller than the symbol 
size. 



the first to flip in the increasing branch of the full hys- 
teresis loop. Therefore, the full hysteresis loop (symmet- 
ric, without EB) will overlap with the minor hysteresis 
loop, as shown in Fig. ||a. In contrast, when the bond 
distribution is distorted, as done in the model presented 
in this work (and which preserves the symmetry proper- 
ties of the Hamiltonian in Eq. [T]), the spins with larger 
ferromagnetic coupling will be the last to reverse in the 
decreasing branch and also the last to be reversed in the 
increasing branch. Thus, minor loops will not coincide 
with the full loops and will easily exhibit large EB as 
indicated in Fig. Ob. 



The justification of the proposed model would require a 
physical explanation for the exchange enhancement phe- 
nomena. Recently it has been suggested that spin-waves 
in the FM/AFM interface could be responsible for an 
enhancement of the exchange coupling p2fl . Here we 
propose a different mechanism based on the existence of 
quenched disorder in the AFM layer. The most common 
case is the existence of antiphase domain boundaries as 
considered in Refs. ||, |2^, [24| and One can assume 
that the exchange interaction between the magnetic mo- 
ments in the FM layer has two contributions: the first 
comming from the direct overlap of the electronic wave 
functions of the atoms in the FM layer and a second aris- 
ing from a superexchange interaction through the overlap 




FIG. 7: Dependence of the relative exchange bias 

{\H e b\) I (AH) on the amount of disorder a for Je = 20 and 
different values of / as indicated by the legend in Fig. ^. Data 
have been obtained by averaging 1000 realizations of a system 
of size L = 50. The lines are a guide to the eye. Statistical 
error bars are smaller than the symbol size. 
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FIG. 8: Hysteresis loop exhibiting exchange bias obtained 
from a numerical simulation of a L — 50 system with Je — 20, 
/ = 0.05 and a = 3. The external field has been cycled 
between H — ±4.5. 



with the electrons in the AFM layer. The existence of 
quenched disorder in the AFM layer can modify this sec- 
ond contribution of the exchange interaction giving rise 
to the exchange enhancement. Since the energies asso- 
ciated with the broken AFM bonds can be higher than 
the FM exchange energies (for instance due to the exis- 
tence of strong anisotropy) , it is plausible to imagine that 
the defects in the AFM can influence the FM exchange 
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FIG. 9: Schematic examples of the hysteresis loops obtained 
for a model with strongly enhanced ±h random fields and for 
strong enhanced random bonds. In both cases the full loop 
is symmetric, but the minor loop only exhibits EB for the 
enhanced random bonds. 



interactions. 

We would like to give a possible mathematical formula- 
tion for such a physical mechanism within the framework 
of lattice models. Let us consider that the interaction be- 
tween two neighbouring spins Si and Sj of the FM layer 
is given by: 



E; 



-JoSiSj 



K SiSjCTi/jj 



(8) 



where at and <jj are the spin variables describing the 
magnetic moments in the AFM layer that sit exactly be- 
low the Si and Sj spins of the FM layer. The constant 
Jo > accounts for the exchange interactions in the free 
FM layer and K > accounts for the coupling between 
the two layers. This coupling term can include, in an 
effective form, the interactions with many atomic layers 
in the AFM. Figure [l(] shows a schematic representation 
of the spins of the FM and the AFM layer in two sit- 
uations: (a) corresponds to the normal situation for an 
ordered AFM layer, whereas (b) corresponds to the case 
in which the AFM exhibits an antiphase domain bound- 
ary. 



FIG. 10: Schematic representation of two spins Si and Sj 
of the FM layer on (a) an ordered AFM layer and (b) on an 
antiphase domain boundary. 



In case (a), the energy of the i,j pair is = 
— (Jo — K)SiSj, whereas in case (b) (or on any antiphase 
boundary) the energy is given by Eij = —(Jo + K)S{Sj. 
The ratio between the two exchange constants is Jg = 
(Jo + K)/(Jo — K), which can be much larger than 1 
when K is close but smaller than Jq. 



V. COMPARISON WITH EXPERIMENT 

The idea that exchange bias could be a minor loop ef- 
fect was already suggested to explain the influence of the 
AFM spin- flop in EB |26| |. Moreover, it has been shown 
experimentally that the uncompensated spins which have 
been suggested to result in exchange bias and which 
could result in the proposed exchange enhancement, can 
be reversed at high enough fields [JIT] (much higher than 
those usually used in exchange bias studies). Hence, in 
this case EB could also be considered as a minor loop 
effect to a certain extent. 

Moreover, the peak observed in the behaviour of the 
coercivity with Je (see fig. |J) could be correlated with 
the experimental observation of peaks in the coercive 
field H c (tAFM) (where tAFM is the thickness of the AFM 
layer) close to the critical AFM thickness for the onset 
of EB or in H C (T) close to the Neel temperature (see 
figures 10 and 13 in Ref. ||). One could consider that 
as the AFM thickness or the temperature increases, the 
FM-AFM coupling, which causes the enhancement of Je, 
decreases and is consequently the origin of the observed 
H c increase, in agreement with theory. Similarly, the 
changes in (Heb) and (Ai?) with could also be cor- 
related with the behaviour experimentally observed in 
Fe/MnF 2 under a cooling field fll^] . 

Another interesting result of this model is the be- 
haviour of EB with the amount of disorder er, as seen 
in Figs. |^ and [?]. Experimentally the role played by dis- 
order in EB is not clear. In some cases there is evidence 
that increasing disorder increases Heb, whereas in other 
cases the opposite effect has been found [|28[ . Our model 
is able to explain both possibilities. Moreover, it is note- 
worthy that experimental non-monotonic dependence of 
Heb on disorder (e.g. roughness, irradiation damage or 
structural disorder) have also been reported ^ ^9|. [w| . 

Finally, another remarkable result from our model is 
that the loops exhibiting EB also exhibit a vertical shift 
in the magnetization axis. This effect has been observed 
experimentally p3, |28|, |3l], 32, 53|. However, since the 



fraction / is very small, the shift can also be very small 
and in some cases, would be difficult to observe experi- 
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mentally. Fig. |Tl| shows the dependence of the magneti- 
zation shift rag on Heb for different values of /. The 
points correspond to different values of a ranging from 
to 2. In view of these results, it would be interesting 
to measure such displacements in different experimental 
systems. 
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by an external field with metastable dynamics. The key 
ingredient is the existence of a tiny fraction of ferromag- 
netic bonds, which are strongly enhanced. This creates a 
long reversible plateaux in the decreasing branch. When 
this pseudo-saturated state is reached, the reversal of the 
field gives rise to an EB loop. Different properties of these 
loops have been computed as a function of the fraction of 
enhanced bonds and the amount of disorder in the sys- 
tem. We have suggested a possible physical mechanism 
to justify the existence of such an exchange enhancement 
in a FM layer on a AFM layer with antiphase domain 
boundaries, based on the existence of superexchange cou- 
pling between the two layers. 

The main conclusions of the paper are: (i) EB is due 
to a minor loop effect; (ii) compensated AFM layers can 
exhibit exchange biased loops with a concomitant mag- 
netization shift; (iii) many experimental phenomena re- 
lated to exchange bias such as: peaks in the coercive 
field, magnetization shifts, marked coercivity increase or 
non-monotonic dependence of Heb on disorder can suc- 
cessfully be reproduced with this model. 
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FIG. 11: Magnetization shift tub versus Heb for different 
values of / as indicated by the legend. The points correspond 
to systems with different values of a, ranging from to 2. The 
lines are a guide to the eye. Statistical error bars are smaller 
than the symbol size. 
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